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Abstract 

The spatial structure of wave functions of anomalously localized states (ALS) 
in disordered conductors is studied in the framework of the cr-model approach. 
These states are responsible for slowly decaying tails of various distribution 
functions. In the quasi-one-dimensional case, properties of ALS governing 
the asymptotic form of the distribution of eigenfunction amplitudes are in- 
vestigated with the use of the transfer matrix method, which yields an exact 
solution to the problem. Comparison of the results with those obtained in the 
saddle-point approximation to the problem shows that the saddle-point con- 
figuration correctly describes the smoothed intensity of an ALS. On this basis, 
the properties of ALS in higher spatial dimensions are considered. We study 
also the ALS responsible for the asymptotic behavior of distribution functions 
of other quantities, such as relaxation time, local and global density of state. 
It is found that the structure of an ALS may be different, depending on the 
specific quantity, for which it constitutes an optimal fluctuation. Relations 
between various procedures of selection of ALS, and between asymptotics of 
corresponding distribution functions, are discussed. 

PACS numbers: 71.55.Jv, 7L20.-b, 05.40.-Fj 
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I. INTRODUCTION 



Mesoscopic fluctuations of various pliysical quantities in disordered systems liave been 
intensively investigated during the last decade 0. These fluctuations originating from the 
quantum coherence of wave functions are typically much stronger than what usual statis- 
tical considerations would predict. In particular, it was found by Altshuler, Kravtsov and 
Lerner that distribution functions of conductance, density of states (DOS), local den- 
sity of states (LDOS), and relaxation times have slowly decaying logarythmically-normal 
(LN) asymptotics at large values of the arguments. These results were obtained within the 
renormalization group (RG) treatment of the field-theoretical a-model describing the low- 
momenta physics of the problem. The validity of this RG approach is restricted to 2D and 
2 + e-dimensional systems, with e <^ 1. 

On the other hand, the conductance, LDOS and relaxation times fluctuations in strictly 
ID disordered chains, where all states are strongly localized, were studied with the use of 
Berezinski and Abrikosov-Ryzhkin techniques IQ-Q. The corresponding distributions were 
found to be of the LN form, too. It was conjectured on the basis of this similarity that 
even in a metallic sample there is a finite probability to find "almost localized" eigenstates. 
Precise meaning of these words remained however obscure. 

More recently, the interest to statistical properties of eigenfunctions in disordered and 
chaotic systems started to grow. On the experimental side, it was motivated by the possibilty 
of fabrication of small systems (quantum dots) with well resolved electron energy levels ^ . 
Conductance fluctuations of such a dot are related to statistical properties of wave function 
amplitudes [§]. Besides, the microwave cavity technique ^ allows to observe experimentally 
spatial fluctuations of electromagnetic wave amplitude in chaotis or disordered cavities |]10 



On the theoretical side, the recent progress is due to the observation that the statis- 
tics of eigenfunctions intensities can be very efficiently studied with making use of the super- 
symmetry technique. This allows one to reformulate the problem in terms of the supermatrix 
cr-model It was found that the zero-dimensional approximation to this cx-model re- 



produces the Random Matrix Theory (RMT) results. In the case of quasi-lD geometry, 
the model was solved exactly by means of the transfer-matrix method |]T^, which allowed 
us to calculate various statistical properties of eigenfunctions |[TT|jr?| -pO| . It was found that 
the distribution function V{u) of eigenfunction intensities u = |^/'^(r)| deviates slightly from 
its RMT form for not too large u, but decreases much slower than RMT predicts, in its 
far "tail". In > 1 dimensions, the results for V{u) were obtained by means of perturba- 
tive method |2y] and saddle-point approximation |21|, and are qualitatively similar to the 
quasi- ID behavior. 

It should be noted that the saddle-point method for the supermatrix a-model used in 
Ref. 1^ and suggested previously by Muzykantskii and Khmelnitskii ||2^ turned out to 
be a powerful tool to study the asymptotic "tails" of various distributions. In particular, 

m 



to find 



it was used in to study the long-time relaxation phenomena, and in 

the asymptotical behavior of the distribution of LDOS. The obtained decrease rates are 
much slower than those given by the perturbation theory. In 2D, the far asymptotics are 
of LN type [^,0, in agreement with the RG results of p|. The asymptotic behavior of 



the distributions was again attributed to the effect of almost localized states p2| , pO| , |2T 
Moreover, it was conjectured in ^^at the form of the saddle-point solution of the 
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cr-model directly describes the electron density of such a state. 

The purpose of this paper is to study in most detail the spatial structure of the anoma- 
lously localized states (ALS). The paper is organized as follows. In section ^ we review 
briefly the previously obtained results for the distribution of eigenf unctions intensity. In 
section [TTTI , which is the central one for the paper, we study the case of a quasi-lD geom- 
etry of the sample. We present the exact calculation of the properties of ALS responsible 
for the asymptotics of V{u). In subsection [111 A 1| the average intensity (|^^(?")|) of such a 
state is found for a sample with the length L ^ C,, where ^ is the localization length. The 
state with an anomalously high local intensity |'?/'^(0)| = u is characterized by an effective 
localization length ^^.j ~ i^/uSY^"^ (here 5* is the wire cross-section), where most of its nor- 
malization is concentrated. Outside this central region, {\ip'^{r)\) decreases as Finally, 
in the vicinity of the observation point there is a sharp change of d'?/'^!) from |-?/'^(0)| = u to 
d"^^!) — ('^/^*S')^/^/2. In subsection pi A 2| we repeat this calculation in the case of metallic 
sample with L <^ ^. We find that ALS have exactly the same form of dip"^ (r) | ) , if -u 3> ^/ SL"^. 
The latter condition means that ^g/ ^ L, i.e. the state plays a role of a localized one even 
from the point of view of a short sample with the length L. 

In subsection piB| we study the fluctuations of intensity |V'^(^)| of an ALS. We find that 
these fluctuations are of the RMT form. The narrow region around the observation point 
where the quasi-jump of occurs, is an exception; here the fluctuations are suppressed. 

Comparing our results with the corresponding saddle-point solution ||2l|, we conclude 
that the latter describes exactly the form of the average density (|'?/'^(r)|), up to a normal- 
ization factor and without the quasi-jump around the observation point r = 0. This allows 
us to generalize the results to the spatial dimension d > 1 (section |V|). In section [V| we 
discuss and compare the spatial structure of ALS corresponding to asymptotical behavior 
of various distribution functions. In Sec.|V5 we discuss briefly some special features of the 
ALS located close to the sample edge, and dependence of the results on the symmetry of 
the ensemble. Our results are summarized in section |V11| . 



II. DISTRIBUTION OF EIGENFUNCTION AMPLITUDES: OVERVIEW OF 

RESULTS. 

In this section, we have collected some results for the eigenfunction intensity distribution 
function obtained with use of the supermatrix a-model formalism. The cx-model action 
reads 



S[Q] = -f^ I d'^r Sti 



TTUqD 



iyQY - TTiyor]AQ 



(2.1) 



Here Q{r) isa4x4or8x8 supermatrix field for the cases of unitary symmetry (bro- 
ken time reversal invariance) and orthogonal symmetry (unbroken time reversal invariance) 
respectively. We will label the formulas for the first case by the index U, and for the 
second case by the index O. The matrix A is defined as A = diagjl, 1,-1,-1} (U) or 
A = diag{l, 1, 1, 1, —1, —1, —1, —1} (O) and the coefficient (3 is equal to /3 = 2 (U) or /5 = 1 
(O). Furthermore, Str stands for the supertrace (trace over bosonic components minus trace 
over fermionic ones), D is the diffusion constant, rj is the level broadening (or, else, (— 2z)x 
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frequency) and is the mean DOS. We do not go into details of the supersymmetric for- 
mahsm here, which can be found e.g. in [^,2^,1^. Let us define now the function Y{Qq) 
as 



Y{Qo) 



Q(ro)=Qo 



DQ{r)exp{-S[Q]} . 



(2.2) 



For the invariance reasons, the function Y{Qo) turns out to be dependent in the unitary 
symmetry case on the two scalar variables 1 < Ai < oo and — 1 < A2 < 1 only, which are 
eigenvalues of the "retarded-retarded" block of the matrix Qq. Moreover, in the limit r] ^ 
(at a fixed value of the system volume) only the dependence on Ai persists: 

Y{Qo) = r(Ai, A2) ^ Ya{2nuovXi) (2.3) 

The distribution V{u) of the eigenf unctions intensities u = |'?/'(roP is then given by pJ]JT5 ] 



V(u) 



Vdu^ 



Ya{u) 



Vdu^ 



u 



2'Ki'Q'q ' 



(2.4) 



In the case of orthogonal symmetry, Y{Qq) = F(Ai,A2,A), where 1 < Ai,A2 < cc and 
— 1 < A < 1. In the limit 77 — 0, the relevant region of values is Ai ^ A2, A, where 

Y{Qo) ^ YaiTTUoVXl) (2.5) 

The distribution of eigenfunction intensities is in this case expressed through the function 



Ya as follows [11 



1 



TtVu^/'^ Ju/2 

ttVu^/"^ du^ Jo z^l 



dz{2z - 

°° dz 



u 



-1/2 



rf2 



Ya{z + u/2) 



YJz) 



(2.6) 



In the diffusive sample, typical configurations of the Q-field are nearly constant in space, 
so that one can approximate the functional integral (|2.2| ) by an integral over a single su- 
permatrix Q. This procedure, which makes the problem effectively zero- dimensional and is 
known as zero-mode approximation, gives 



and consequently. 



Ya{z) 
V{U): 



-Vz 



Ve 



-uV 



2ttu 



{0,U) , 



(2.7) 

(2.8) 
(2.9) 



which are just the RMT results for the Gaussian Unitary Ensemble (GUE) and Gaussian 
Orthogonal Ensemble (GOE) respectively . 

For the case of quasi-lD geometry one can solve the problem exactly, by evaluating 
the integral in ( p.2| ) with making use of the transfer-matrix method ||lT|. As a result, one 
gets an expression for V{u) depending on a scaling parameter x = L/^, where L is the 
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sample length, ^ = 27ri'oSD is the localization length, and S is the transverse cross-section 
|TT| , [T9| , ^ . In particular, in the case of a short (metallic) sample, x <C 1, we find for the 
distribution of normalized intensity y = uV: 



-h/2 



l + ^{2-Ay + y'') + ... 
1 + f (|-3y + ^)+.. 
-y + fy^x + ...]} ; 



P(y) ~exp{f 

V{y) ~ exp \-2l3^JyJx 



y < 

y < X-'/' 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



Here the coeficient a is equal to a = 2[1 — 3L_L_(_/L^], where L is the sample length and L_, 
are the distances from the observation point Vq to the sample edges. If one is interested 
in the distribution function V{y) averaged over the position of the point Vq within the 
sample, one has to replace a in eqs.( |2.1UD , ( |2.11D by its average value equal to unity, and 
in eq.( p.l2| ) by its maximum value equal to 2. Such averaged distribution functions were 
considered in our previous publications, Ref. ||TT1JT9| , |20[| . It is straightforward, however, to 
generalize the derivation and to get the formulas ( p.lO| )-( p.l2D for the position-dependent 
distribution functions. Note that g = 1/x is for small x just the dimensionless conductance 
of the sample: g = 2nh'QDS/L = G/i^e^/h). We see, that for not too large amplitudes y, 
Eqs. (|2.1CI| ), (|2.11|) are valid which are just the RMT results with relatively small corrections. 
In the intermediate region ( |2.12 ), the correction in the exponent is small compared to the 
leading term but much larger than unity, so that Viy) ^ VuMriy) though lnP(?/) ~ 
InVjiMTiy)- Finally, in the large amplitude region, ( |2.13D , the distribution function V{y) 
differs completely from the RMT prediction. We will present eq. ( |2.13| ) in a more precise 
form in Sec.|l|. Note that it is not valid when the observation point is located close to the 
sample boundary, as will be explained in Sec. |VI A| . The asymptotic behavior (12.131) is the 



same as in the case of insulating sample, L ^ ^, when the distribution V{u) takes the form 



L 



Kti2JuSO+m2JuSO 



2i'^S Ki{2^/^) 



(2.14) 
(2.15) 



This allowed us to conjecture in Ref. that the asymptotic behavior ( [^.13| ) is controlled 
by the probability to have a quasilocalized eigenstate with an effective spatial extent much 
less than ^. We will prove this rigorously in the present paper. 

In d > 1 dimensions the problem can not be solved exactly anymore, and some ap- 
proximate methods of evaluation of the integral ( p.2|) are necessary. For moderately large 
amplitudes y ^ a perturbative treatment of non-zero modes of the a-model in the 



weak localization region is possible ||20|. Here k is the usual parameter of the perturbation 
theory M defined as 



E 



1 



2-nvQDq^V 



(2.16) 



the summation goes over the eigenmodes of the Laplace operator in the sample. The result 
for Viy) is: 
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e-y[l + K{2-Ay + y^) + ..] 



1 P ^ ^y' 



+ 



(2.17) 
(2.18) 



In particular, in quasi-lD systems k = l/Qg, and Eqs. (|2.1CI|) , (|2.11| ) are reproduced. For 
d > 2 the corresponding sum over momenta q diverges at large \q\ and is to be cut off at 
l^^. This gives k = ^^2^^j~, ln(L//) for d = 2 and k ~ {kplY d = where kp is the 



Fermi momentum. 

In the region of large amplitudes y ^ k^^/^, Eq.( |2.17D , ( |2.18| ) loose their validity. In this 
region, the distribution P(|/) can be found |2ll] by using the saddle-point approximation. 
For (i = 2, 3 the result is 



P(l/) ~exp{f(-y + /€y2 + ...)| 
P(y)~exp{-£ln'^(/s:y)} , 



K ^ ^ y ^ K, 
y ~ 



(2.19) 
(2.20) 



In the next two sections, we will study in detail the structure of the states, which are 
responsible for the far asymptotics ( p.l3| ), ( |2.20 ). For definitness, we will consider the unitary 
symmetry case; generalization to the orthogonal ensemble is completely straightforward, and 
we will just quote the (minor) changes in the results in Sec. [VI B| . 



III. ANOMALOUSLY LOCALIZED STATES IN QUASI ID SAMPLES 

A. Average form of the eigenfunction corresponding to an anomalously high local 

amplitude. 

In this section, we study the average intensity {\'4^'^{f)\)u of a wave function with the 
condition that |'?/'^(0)| = u. For convenience, we have put the observation point in the 
coordinate origin; the sample edges are at r = — L_ and r = L_|_; the sample length is 
L = L_ + L+. Formally, {\'ip'^{r)\)u is defined as 

(IV^V)I). = ^|^, (3.1) 

where 

V{u) = -^{J:S{\M0)\' - u)6{E - E^)) (3.2) 

Q{u,r) = \iY.\i'^i'')\'^i\^M?-n)5{E~E^)) , (3.3) 

where V^q, are eigenfunctions of the Hamiltonian, and E^ are the corresponding eigenvalues. 
The here defined function V{u) is just the distribution function of local intensity discussed in 
the preceding section. As has been mentioned above, it was calculated in the quasi-lD case, 
with making use of the a-model representation, eqs.(|2.4]), (f2.6|), and the transfer-matrix 



method [11,191. The result reads 
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V{u) = l|l{w^«(«5e,r+)iy«(«5e,r_)} , (3.4) 
where the function W^^\z^t) satisfies the equation 

and the boundary condition 

W^^\z,Q) = l (3.6) 
The solution to eqs. (|3.5|) , ( p.6|) can be found in terms of the Lebedev-Kontorovich expansion 

W^^\z, t) = 2^1/2 \k,{2z^'^) + - r du-^ sinh ^KU2z'/')e-"-^A (3.7) 

I TT JO 1 + Z/^ 2 J 



Now we turn to evaluation of the function Q{u,r) defined in eq.( p.3| ). Detailed exposition 
of the method used can be found in Ref. [|19|, and we will mainly follow the notations of 
this paper. We start from expressing the moments (|'?/'(r)p|'?/'(0)p'^) in terms of the Green's 
functions 



where 



hm— — (27^)^(G?j(0,0)GA(r,r)) , (3.8) 



GRAri.r2) = {rr\{E-H±ir^)-'\r2) (3.9) 



Here H is the Hamiltonian which consists of the free part Hq and the disorder potential 
f/(r): 

H = Ho + U{r)- H, = -Lp^- (3.10) 

2m 

the latter being defined by the correlator 

{U{r)U{r')) = -^6{r-r') (3.11) 

Next, we write the product of the Green's functions in terms of the integral over a supervector 
field $ = {Si,S2,Xi,X2)- 



Gl{0,0)GA{r,r) = ^ / /^$t(5i(0)S*(0))'^S2(r)S*(r) 

q. J 

X exp |z J dr'<l>\r')A^/\E + ir]A - H)A^/^<l>{r')^ (3.12) 



The following steps are [IS 
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i) averaging over the disorder; 



ii) introducing a 4 x 4 supermatrix variable R^y{r) having the symmetry of the tensor 
product <l>^(r)$J,(r) ; 

iii) integrating out the $ fields ; 

iy) using the saddle-point approximation which leads to the following equation for R: 

^ ' ^ (3.13) 

(3.14) 



g{ri,r2) = {r^\{E-H^-R)-^\r2) 



The relevant set of the solutions (the saddle-point manifold) has the form: 



R = a-I Q 



(3.15) 



where J is the unity matrix, a is certain constant, and Q belongs to the coset space t/(l, 1|2) 
and satisfies the condition = 1. Finally, we find 



(i^wn^(o)r^) 



— \Ycv.{2T^v,ri)^ \ DQ Q?i,,,(0)Q22,56(r) 



+ ?^^^^Q?li(0)^12,66(0, r)^21,feb(T', 0) 



(3.16) 



where S\ff\ is the cr-model action defined in eq. 

Taking into account that Qiv^ varies slowly on a scale of the mean free path /, we have 
for Iri — <^ / 



(3.17) 



where 



(3.18) 



is the average one-particle Green's function. In the opposite limiting case \r\ — ^ /, the 
Green's function girx^r-i) vanishes exponentially. Thus, we get 



2V 



22,bb\l )^ 



-S{Q} 



r:^l (3.19) 



mr)?\m?') ^ -^rM'^nmY / DQ g?i,,,(o)Q22,66(o) 



+ gA;,(r)g?-i(o)gi2,b6(o)g2i,6fe(o) 



-S{Q} 



r < I 



(3.20) 



where 
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(3.21) 



depends on short-scale dimensionality d of the sample [^. In particular, for a strip (d = 2) 
or a wire {d = 3) we have 



2 ' 



(3.22) 
(3.23) 



respectively. 

Now we evaluate the Q-integrals in eqs. (|3.19|) , ( p.20| ) using the transfer-matrix method 
I9| . We first consider the r ^ / case. We will assume that the transverse size of the 



wire is much less than the effective localization length of the ALS, which we will find to be 
^e/ ~ This means that we are indeed in the quasi-lD regime, when considering the 

ALS structure. We assume for definitness that the point r is located to the positive direction 
from the coordinate origin: —L_ < < r < L^. We get then: 



dzz'^-^W^^^z, T_)W^^\z, Ti, Ta) 



(3.24) 



where r_ = L-/^, ti = r/^, T2 = (L+ — r)/^; the function W'^^^z, r) is defined by eqs.(p3| 
(^.7|), and the function W^'^\z,ti,T2) satisfies the same equation 



Oti \ dz^ 

and the boundary condition 

The solution of eqs.(|3:25|), ( CT ) is H 

W^''\z,T^,T2) = 2v^ r dph{u,T2)Ki,i2^z)e-'^ 

Jo 



(3.25) 



(3.26) 



V sinh(7rz/) 

2^i^2 



dtK,,(t)W^^\t^l^, 



^2) 



(3.27) 



Substituting here the formula ( |3.7|) for W^^\z^T2) and evaluating the integral over z, we 
reduce eq. (|3.27|) for 6(1/, T2) to the form 



K^, ^2) 



r 



l + iv 



+ 



V sinh(7rz/) 

16^^2 

V sinh(7rz/) f dvi Ui 



27r3 



sinh 



(l + u' 



r 1 + 2 



1^ + 1^1 



r 1 



.z/ — z/i 



e — (3.28) 



Equations ( |3.24| ), (|3?^ ), (|3.27|) , (|3.28|) constitute the final result for the moments 
(|-?/'(r)p|?/'(0)p'') at r ^ /. Now we can restore the function Q{u, r): 
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Q{u,r)^{5{\m?-nmr)\') 



d 



ViSdu 



u 



r > / 



In the opposite case, r <^ /, we find from eq.( |3.2CI| ) 



d 



d 



QM = - Ur) \u— + -]--) YM 



V 



du^ du 



du 



(3.29) 



(3.30) 



wliere tlie function Ya{u) was defined in eq.( p.3|) . This formula is vahd for any sample, which 
is locally c?-dimensional. In the case of the quasi-lD geometry we get 



1 

V 



V dw^ du 



du 



r<^l (3.31) 



Substituting these results along with the formula ( |3.4| ) for V{u) in eq.(p.l|), we can find 
the average density {\ip'^{r)\)u- It is possible to check that it satisfies the normalization 
requirment / d'?/'^ (*')!)« = 1, as it should be. In the following two subsections, we will study 
the form of (|V'^(^')|)u for an ALS in the insulating (L ^ ^) and the metallic {L ^ ^) limits. 



1. Insulating sample (L ^ 

In this subsection, we analyze the above general results for {\ip'^{r)\)u in the limit of an 
infinite sample length. More precisely, we will assume that the distances from the 

point r = to the both edges of the sample are much larger than ^, so that we can consider 
the limit t_, r+ oo. We have then 

W^^\z, oo) = 2v^iri(2v^) (3.32) 

and the distribution function V{u) takes the form ( p.l4| ). Furthermore, eqs.( p.27| ), (p.28|) 
reduce to 

W^''\z,n, oo) = I4~z /"dz/^^^^^^tanh^iri,(2v^)e-^"i (3.33) 

JO 8 2 

Typical values for the intensity m of a localized state are u % ^/SC,- We are interested 
however in anomalously high u ^ l/5'^. It is seen from eqs. (|3.4| ), (|3.7|) , (|3.29|) , (|3.33|) , that 
this corresponds to large values t = 2^JuS^ ^ 1 of the argument of the modified Bessel 



functions Ki{t) and Kiy{t). The corresponding asymptotic formulae read 



^iW^/fe-*, t»l (3.34) 
ir,,(t)~yjexp|-t-t/(^^)| , Kt, v<t (3.35) 



f{w) = y/l — + w arcsin w — 1 = — + 0{w'^) 
The asymptotics of the distribution function V{u) has the form 



2 
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VV{u) ~ 47r5=^/2^3/2^-i/2g-4y^ 



(3.36) 



As to the function r), we consider separately three regions of the distance r: 

a) r 3> ^, i.e. ri = r/^ > 1. 
The integral in ( |3.33|) is then determined by the region ~ ri yielding 



W^'^\z = tV4,ri,cx)) 



8^2 



(3.37) 



where we have denoted t = 2yi. Consequently 



VQ{u,r) 



4m 



3/2 -4yM5C-ri/4 



(3.38) 



b) / < r < ^. 
We have 



1 irf 

W^^\z = 1" / A, Ti, oo) c^-\ —e-'-^^ I ^ / rfz/z/(l + z/2)tanh 

8 V 2 JO 



TTi/ r 2 1 ri 



so that 



" f5/2 

2 (rit + 2)^ 



l + rivW)2 



(3.39) 



(3.40) 



c) r -C /. 

In this region an additional term proportional to the Friedel function kd{r) appears in the 
expression for Q(m, r), see eq.( p.31|) . We get 



VQ{u,r) = 2-KiSe 



1 + 2JuSikd{r) 



(3.41) 



Since fcrf(r) ~ 1/ {kprY'^ for r 3> k'^^^ the second term in square brackets is much larger 
than the first one at r <C tq ~ ^^^(mS'^)-^/'^^'^"^)] and is negligible at r ^ tq. 

Substituting now eqs. (|3.36| ), ( |3.38|) , ( p.40|) , and (|3.41| ) in eq. (|3.1| ), we find the following 
spatial structure of the ALS with |'?/'^(0)| = u: 



{\^\r)\)^ = ^u~i/25"3/2r-3/2e-/4? 



(I^^(r)l). 



{mr)\)u 



u 



1/2 



2 \iS 

1 /«\l/2 



2 ' 



2 VS-S 



1 + 2^^kd{y 



/ < r < ^ 
r < / 



(3.42) 
(3.43) 

(3.44) 



We see from eqs.( p.4^ ), ( |3.43| ), ( p.44| ) that the eigenfunction normalization is dominated 
by the region r ~ ^e/, where ^e/ ~ y^^/'uS ^ ,^ plays the role of an effective localization 
length. In the region S,ef ^ r ^ ^ the wave intensity falls down as 1/r^, and transits to the 
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conventional localization behavior at r ^ ^. Therefore, the appearance of an anomalously 
high amplitude |'?/'^(0)| = u ^ ^/S^ is not just a local fluctuation, but rather a kind of 
a cooperative phenomenon corresponding to existence of a whole region r < C,ef with an 

unusually large amplitude IV'^I'^)! = i\AV6S' ~ 

Let us emphasize once more that what we have calculated is the average value {\ip'^{r)\) 
of the eigenfunction intensity with the condition |V'^(0)| = u. It is natural to ask now what 
are its fluctuations. This question is addressed in section |III B| . We will also explain there 
what is the reason for the "quasi-jump" of (I'lp^l) from |'?/'^(0)| = m to {\'ip'^{r)\) = ^^u/C^S 
at r ~ /. 



2. Metallic state (L <^^). 



We will assume that the observation point r = is located somewhere in the bulk of 
the sample, so that both r_ = L^/S, and r+ = are of the same order of magnitude as 

r_,r+~L/e = l/(7<l 



The distribution of eigenfunction intensities is given by eqs.f p^) , ( p. 71) , and its behavior in 
various ranges of the variable y = uV is indicated in eqs.( |2.1Cl| )-( |2.13|) . We will study the 
structure of the ALS responsible for the far asymptotics (|2.13|) in the region uV 3> g. As 
is seen from eqs.( p.4| ), ( |3.7| ), this corresponds to a large value of the argument t = 2z^l'^ = 
2(uS^Y^'^ ^ g of the modified Bessel function in eq. (|3.71 ). Under this condition, the integral 



in eq.( |3.7| ) can be evaluated via the stationary point method with use of the asymptotic 
expression ( p. 35 ) for the modified Bessel function: 



= iy2Fexp{-t + =!-t/(^)) . 



t > 1/r > 1; 



(3.45) 



w 



+ ... , 



w < 1 



The distribution function V{u) has therefore the form 



V{u) = ^J^ 



exp <^ -4a/<5 + — ^ + 



4L, 4L_ 



f 



i|^^exp{-4;;^^+£( 





+ 



(3.46) 
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This is just the formula (|2.13| ), but written now with full accuracy with respect to the 
subleading factors. 

To calculate Q{u,r), eq.( p.29| ), we have first to evaluate the function W^'^\z,Ti,T2), 
eq. (|3.27|) . The contribution to it from the first term in eq. (|3.28| ) was calculated in the 
preceding subsection, where the insulating regime was considered. We will find that in the 
metallic regime the second term in eq. ( |3.28|) gives a much larger contribution, so that the 
first one can be neglected. Substituting the second term of eq.( |3.27| ) into ( |3.28D , we find 
after simple algebraic transformations: 



t 

Air 



smh 



z/sinh(7rz/) 



Z/2 - Z/2 



2 cosh(7rz/i) — cosh(7rz/) 



(3.47) 



Analysis of the double integral in eq. ( |3.47D shows that one should distinguish between the 
two possible siuations: 
i) Tl + f < r| ~ l/g\ 

This corresponds to a very large amplitude ^ ~ |^ ^ V^^g^. The leading contribution 
to the integral in p.47| ) comes from the region u ^ Ui. To check this, let us assume that 
u ^ Ui and collect the exponential factors in eq. ( 13.47] ): 



SO that the characteristic values of the variables are z/ ~ (ri + 2/i(:)"^/^ and vi ~ 1/t2, con- 
firming the consistency of our assumption in the considered range of parameters. Collecting 
all the prefactors, we get 



2ti\ 2 (l/t + ri/2)2 



^7r2/4T+-t . 



ri + T2 



(3.48) 



Note that we have omitted the corrections of the type tf{TT/tT^) in the exponent, since they 
are small in the considered case g"^. Substituting eq. ( |3.48|) in eq. (|3.29|) , we find 



8 ^fL^7- 



(l/t + r/20- 



t = 2JuSi , (3.49) 



and finally 



mr)\)^ 



1 (u^ 



1/2 



(3.50) 



ii) < Tl + f . 

In this case the integrals in (|3.47|) are dominated by the domain z/i ~ z/ ^ 1. Introducing 
the new variable z/_ = z/ — z/i , we reduce eq. ( |3.47D to the form 



13 



1 t 



8 V 27r 
X exp 

Evaluating the integrals, we find 



dvdv^ 



t-tf 



sinh(7rz/_/2) 



TTZ/ 

T 



-n 



-7-2 



1 /2t 
Att]] tI 



C 2, 1 



2r+ ; 
X exp 

where ({p, z) is the generalized Riemann's zeta- function: 

oo 

ap,z) = Y.{z+k)-^ 



TX 



2 



4n 



2r^ 



(3.51) 



(3.52) 



fc=0 



Substituting this in eq.( p.29|) , we get 



Q{u) 



C 2, 1 



X exp 



2r+ . 

+ 



C 2, 



2t-t/ 



2t^ 



71 



tf 



^4r+ 4r_ 

Therefore, the average local intensity of the ALS is given by 

1 



7r 



t = 2duSi (3.53) 



(l^'WI)n 



85L 




C 2, 



2L, 



C 2, 1- 



2Li 



(3.54) 



In particular, in its "core", r <^ L+, the ALS intensity has again the form ( [3.501 ). 

Comparing eqs.( |3.50| ), ( p.54| ) with the result (|3.43| ), we conclude that in its central "bump" 
the ALS in the metallic sample has exactly the same spatial structure as in a long (insulating) 
one. The condition Vu ^ g, under which the asymptotical behavior (|2.13|) , (|3.46|) is valid, 
acquires now a very transparent meaning. This is just the condition that the effective 
localization length of an ALS, ^^f = is much less than the sample size L. Indeed, 

ief/L = ^^/uSL"^ = ^g/uV. In the "tail", r L ^ ^e/, the form of the ALS intensity is 
slightly modified by the boundary of the sample, see eq.( p.54|) . 

All the above calculations in this subsection were valid for the domain r > I. As to 
the region r < /, we can easily satisfy ourselves using eqs.( p.3lD and ( p.45| ) that the result 
( p.44| ) holds. Therefore, the "quasi-jump" of {\ip'^{'r)\)u at r <^ / has the same form, as in 
the insulating regime. 



B. Fluctuations of the eigenfunction of an anomalously localized state. 

In the section [III A| we have calculated the average intensity {\ip'^{'r)\)u of a quantum state 
which has an anomalously high amplitude |V'^(0)| = m in a certain point r = 0. However, this 
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average value does not give full information about the ALS. In the present subsection, we will 
study the fluctuations of the intensity of an ALS fixed by the condition |'i/'^(0)| = u. 

Similarly to the subsection [111 A| (see eq.( p.47| )), we express the moments {\4'{r)\'^P\ip{0)\'^'^) 
in terms of the Green's functions: 

{Y^\Mr)\''\MO)\"'S{E ~ E^)) 

a 

jQ~p (a — \ \](t) — 1 V 

W iJ-lP ^^•.(2r^)«+p-i(Gyo,0)G^(r,r)) , 



lim 



'7-o27rz/oV (g + p-2)! 
The analogue of eq. ( |3.12| ) reads 

^ D<^D^\SMSm)\S2{r)Sl{r)Y 



(3.55) 



p\q\ 



X exp |z J dr'^\r')A^/\E + ir]A - H)A^''^^r') 
Proceeding further as in the section [111 A|, we get for r S> / 



(3.56) 



mr)?''\m?') 



^hm(2..o.)--^^^^^^^^ 



2V v^o' 



(3.57) 



(g + p-2) 

For the quasi-lD geometry the integral can be again evaluated via the method of ||TU 
yielding 



1 



q\p\ 



V{^S)'}+P-^ {q + p-2) 



roo 

{z,Ti,T2) , (3.58) 



where Wj;'^'^ satisfies the same equation ( p.25[ ) as the function ly^^^ = Wi^\ and the boundary 
condition generalizing eq. (p.26[) : 



Wj,^\z,0,r2)=zPW^'\z,T2) 



Defining similarly to eq. ([3.3[) , 

1 



r) = -^(E \Mr)\'''S{\MO)\' - u)5{E - K)) 



we get 



(3.59) 



(3.60) 



— 1-- 



W^^\vuiS, n, T2)W^'\vu^S, r_) 



For the second moment, p = 2, this formula can be simplified 



p> 1 
(3.61) 
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Q2{u,r) 



V{uSO' 

The solution to eqs. (|3.25| ), ( p.59| ) has the form 



2iJ sinh(7rz/) dt ft 



(3.62) 

(3.63) 
(3.64) 



As we have seen above, an ALS in a metalhc sample has essentially the same spatial 
structure as in an insulating one. Physical reason for this is very simple: an ALS with 
the effective localization length C,ef may exist in a sample with the length L ^ C,ef without 
essential modifications. This reason is equally valid for the fluctuations. Therefore, we can 
limit ourselves by studying the fluctuations in the technically simpler case of an infinitely 
long sample. In this case, the function W^^^ is given by eq. (|3.32|) , and eq.( p.64| ) reduces to 



Ki^^ ^2 = 00) 



u sinh(7rz/) 
47r2(2p- 1) 
z/tanh(7ri//2) 
24p-i(2p- 1)1 



T{p + 



l + iu 



Tip 



z/2 + (2p-l)2 u^ + {2p-?>f 



(3.65) 



The integral in eq. (|3.63|) is determined by z/ ^ 1. Therefore, for not too high moments p 
the inequality p <^ v is satisfied; the corresponding restriction on p will be found explicitly 
below. We have then 



6p(i/, 00) 



1 



{2p-l)\ \2 



V 



4p-l 



(3.66) 



and consequently from eq.( p.63|) . 

Substituting this in eq. ( p.61|) , we get 

47rp 



2 (ri + 2/t) 



V Jl \ V J 



2p 



(3.67) 



vu^S 



2p 



(3.68) 



For not too large p <^ the integral is determined by the region f — 1 ^ 1, because of 

the factor exp{—4:y/vu^S}. It can be then estimated as 



Qpiu,r) 



Anp 



{SO 



3/2-p^-l/2 



2p 



I J 



^pliSifl^u-^l^e-"^' 



2(l + riv^)2 



(3.69) 
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Finally, the p-th moment of the ALS intensity is given by 

Qp{u,r) ,\ \RiS 



V{u) 



pi 



2{1 + rJuS/0' 



(3.70) 



In course of the derivation we assumed that p <^ u for typical values of u in the integral 
( p.63| ). These values are 



P 



Tl 



SO that the condition reads 



p < 



1 + ^ 



Pn 



(3.71) 



uSS, 



An ALS is defined by the condition uS^ 1, so that for r <^ ^ we have Pmax ^ 1- 
Comparing eq.( |3.70D with eq.( |3.43| ), we find that in this region 



(|^2^(r)|)„^p! \{\r{r)\)^ 



I < r <^ ^, p <^Pr, 



(3.72) 



This means that the fluctuations of the eigenfunction intensity with respect to its average 
value are of the type 



2(1 + r^uS/ if 
where |$^(r)| is distributed according to the GUE law 



I <r 



(3.73) 



(3.74) 



The approximate result ( p.74| ) holds for not too large |$^| <^ Pmax- 

Now we consider the fluctuations in the vicinity of the observation point, r <^ I. Pro- 
ceeding as in the section [111 A] , we find from eqs.( |3.55D , (|3.56|) 



(|^(r)|^^|^(0)r^) 



2q\ 



P+1- 



p\ 



q\p\ 



jKp-jV-JKq-jV- '{q + p-2)\ 



X dzz'>-^P-^W^^\z,T+)W^^\z,T_) 



(3.75) 



This allows one to restore the joint distribution function of |'?/'(0)p and |'?/'(r)p (see Ap- 
pendix), which looks however rather cumbersome. For this reason, let us return to the 
expression for the moments and consider p = 2. We find then 



viso 



g+1 



l + 2qUr) + ^^^-^kl{r) 



dzz''W'^^\z,T+)W^^\z,T_) 

(3.76) 
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Restoring the function Q2{u,r) defined in eq. (p.60|) , we get 



du I dw^ 



Substituting here (p.32|), we find 



{W^^\uSi,T^)W^^\uSi,T.)} (3.77) 



^ kl{r)u' + 2u^-^fc,(r)(l - k,{r)) + ^(1 - k,{r)f (3.78) 

Therefore, the variance of \^'^\ is equal to 

var„(|^2|)^(|V^^(r)|)„-(|^2(r)|)^ 

4% , r > ro , 

with the scale rg as defined after eq. (|3.41|) . We find that for r <^ tq the fiuctuations are 
suppressed: 

varufl-j/'^fr)!) 1 — kd(r) 1 

' ^ \2 ^ , / « 1 3.80 

On the other hand, for r ^ tq we have the GUE result 

var„(|^2(^)|)^^|^2^^)|^2 ^3^g^) 

More generally, it is possible to check that for r ^ ro the terms containing kd{r) become 
negligible in higher moments of |'?/'('")^| as well, and the GUE-like fiuctuations ( |3.72| )-( p.74| ) 
take place. 

Let us summarize the results of this subsection. We have found that the ALS intensity 
|'?/'^(r)| exhibits in the region tq <C r -C ^ the GUE type fiuctuations ( |3.72| )-( |3.74| ) with 



respect to its average value {\4''^{r)\)u- These fiuctuations are completely analogous to those 
for an ordinary delocalized state in a metallic sample, see eq.( |2.^ ). The difference is that 
in the present case the average intensity {\ip'^{'r)\)u is not uniform in the coordinate space. 
In the region r <^ rg the fiuctuations are suppressed: var„(|-?/'^(r)|) <ti {\'ip^{r)\)l^. This is 
similar to what has been found in pO| for the spatial structure of an "ordinary" delocalized 



state with a moderately large local intensity u = |'?/'^(0)| (when the zero-dimensional formula 
( p.8|) holds and an ALS is not formed). This means that the intensity of a typical ALS is 
in this region close to the average value {\ilj'^{r)\)u, which exhibits the sharp decrease from 
|'?/'^(0)| = u to {\ip'^{r)\)u = at r tq. It is not difficult to understand that this quasi- 

jump has the same sourse as the GUE-like fiuctuations at r 3> tq. One can ask, of course, 
why this short-scale fiuctuation happens exactly in the center of the smooth ALS "bump" 
with a probability close to unity. The answer is as follows. We are studying the states with 
an anomalously large local intensity u, which is an exponentially rare event. There are two 
sources which may favor the formation of such a high intensity: i) formation of an ALS with 
a spatially non-uniform smooth envelope, and ii) short-scale GUE-like fiuctuations. Both 
these mechanisms have exponentially small probabilities to produce an enhancement of the 
intensity by a large factor. The found configuration of {\ip^{'r)\)u (short-scale quasi-jump 
( p.44| ) on top of the smooth configuration ( |3.43| )) represents just the optimal combination 



of the two mechanisms. 
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IV. GENERALIZATION TO HIGHER SPATIAL DIMENSIONS. 



The calculation in the preceding section relies essentially on the quasi- ID structure of 
the sample, so it is not applicable for a sample of higher dimensionality. In this case, the 
asymptotical behavior of the distribution function V{u) was studied by Fal'ko and Efetov 
pT| . These authors used eqs.( p.4|) , (pl6|) and applied the saddle-point method suggested for 
the super symmetric cr-model by Muzykantskii and Khmelnitskii . The advantage of this 
method is that it is applicable for any spatial dimension d. It was conjectured in [0,^ that 
the saddle point solution mimics the spatial form of the ALS. Having at our disposal the 
exact solution for the quasi- ID case, we can check the accuracy of the saddle-point method 
and of this conjecture. Comparing the form of V{u) found in the quasi-lD case by the 
saddle-point method ||2T| with the result of the exact solution ( ^.46[ ), we find a very good 



agreement between them. Furthermore, the saddle-point configuration in the quasi-lD case 
has the form pT,P^| 



AiW _ ^,(.)_,(o) ^ ^ 1 . o<r«L+ (4.1) 



Ai(0) 



1 + r 



Here Ai = e /2 is the eigenvalue of the Q-matrix which has been introduced after eq.( 
Comparing eq.( [4.1| ) with eqs. (|3.43|) , ( |3.50|) , we see that the saddle-point solution nicely re- 
produces the average intensity of the ALS, for r > Z, up to an overall normalization 
factor. Being encouraged by this agreement, we will now use the results of the saddle-point 
study to describe the structure of the ALS in d = 2, 3. 



A. 2D geometry. 

For a 2D disk-shaped sample of a radius L with the high amplitude point r = in the 
center of the disk, the saddle-point solution was found to have the form 

where the exponent < /i < 1 depends on u and satisfies the equation 

^ (4.3) 

We are interested in the asymptotic region uLF' ^ tiv^D ln~^(L/Z), where the distribution of 
the eigenfunction intensity is given by eq. ( |2.2(]| ), and an ALS is formed. Then the exponent 



can be approximated as 

In 

" 21n(L// 



In f In f ) 

^ ~ '-^ (4.4) 
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The lower cut-off scale /* appears in eg. ( |4.2| ) because of the restriction of the diffusion 
approximation on the momenta q of the a- model field: q < Correspondingly, it is 

determined by the condition p2|-p^ 



dr 

which yields ^ fil, so that l^, differs from / by the logarithmic factors (|4.4| ) only. 

Normalizing properly the expression ( [4.2[ ), we find that the average ALS density for 
r > Z^, is equal to 

<l^'(-'l>- - 4^ (9'"{' - si:^%...D ilTT ■ 

The saddle-point calculation of Ref. |2l| assumes that d{r) is constant for r < so that 
eg. ([4 .51 ) gives {\ip'^{'r)\)u — ^^a"^^)^ in this region. However, for very small r < I the average 
intensity {\4''^{r)\)u changes sharply, as we have seen in the quasi-lD case. Indeed, in this 
domain the function Q{u,r) is given by eq. (|3.30|) , where Ya{u) has the following behavior 

El 



Using now eg. ( |3.30| ) for Q{u, r) and eq.( |2.4| ) for V{u) , we get 

(|^2(r)|)„ ^ = [1 - hir) + Aiu)h{rm^'{r = Q\U , r<h (4.7) 

Here the height of the quasi-jump is given by 

Aiu) . ^u± ... y.(.) . g^g .„ .„ . 4,V„Z,, (4.8) 



Combining eqs. (^^ ), (|4.7| ), and (|4.^ ), we get 

{\i^\r)\)^ = —-[1 - h{r) + h{r)A{u)] , r<h (4.9) 

A[U) 

In particular, at r = we find (|^/'^(r)|)u = u, that shows the perfect consistency of the 
whole procedure. 



B. 3D geometry. 



Now, we consider the 3D case. The saddle-point configuration has the form [21,24 



e{r)-9{0) = c(--l] , r>h, (4.10) 



r 
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where the coefficient C is defined by the condition 



Ce 



c 



uV 



(4.11) 



and /* = CI is the cut-off length scale which has the same origin as in 2D. We are interested 
in the asymptotic region uV ^ vqDI, where eq. (|4.11|) yields 



uV 

lHuV/g{l)] , 



(4.12) 
(4.13) 



where g{l) ~ {kpl)'^ is the conductance on the scale of order of the mean free path /. The 
optimal configuration( |irTO| ) has the form 



,0(r)-e(O) 



uV 



^g{l) \ri'[uV I g{l)] 



(4.14) 



This configuration determines the asymptotical behavior of the functions Ya{u), V{n) at 
Vu > g{l): 



V{u) ~ Ya{u) ~ exp < —const g{l) W 



Vu 
W). 



(4.15) 



Note that in order to fix the numerical coefficient in eq.( [4.15|) , one has to go beyond the 
long-wave-length a-model approximation. Such a "ballistic" generalization of the cr-model 



was recently suggested in Ref. |Q. The height of the quasi-jump A{u) is found from eq. (^78| 
to be 



Aiu)r^gil)\n 



2 Vu 



gil) 



(4.16) 



Thus, the density of ALS is 

0(r)-e(O) 

mr)\)u- 



A{u) 



1 

V 



uV 



U/r 



r >U 



(4.17) 



The quasi-jump at small r is given again by eq.( ^.9D , with k2{r) replaced by k^{r). 

The value of (|V'^(^)l)u given by eq. ( |4.17|) exceeds considerably the average density of a 
delocalized state, (|'?/'^|) = l/V at r ^ where 



,1 uV 2 uV 
m — — ~ I m 



9(1) 



9{iy 



For larger distances, r > we have {\ip'^{r)\)u — ^/V. Therefore, in contrast to quasi-lD 
and 2D systems, where an ALS formation is a redistribution of the eigenfunction intensity 
in the whole sample, in 3D it constitutes just a "bump" with the extent of order of on 
top of the usual background density {\ip'^\) = l/V. 
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C. Fluctuations. 



Let us now discuss briefly the fluctuations of the ALS intensity. At small r <^ Z, we find 
in full analogy with eqs.(^:7HD-(F^): 



{\ij\r)\)^^u' 
and consequently, 



1 - Ur) 



(4.18) 



kd{r)A[u) (4.19) 



{mr)\)i 



1 , r > ro 



Here tq is the characteristic spatial scale of the quasi-jump, determined by the condition 
kd{r)A{u) ~ 1. In 2D, using eqs. QO^) and (gj), we find 

In HL/l)) 

In 3D the analogous calculation would give tq ~ l\n.\Vu/ g{iy\ > I. This means that in fact 
tq ~ I, because of the exponential decrease of kd{r) ~ e~^^'' at r > /, which was not taken 
into account in eqs.( p.22D , ( p.23| ). We see from eq.( |4.19| ) that at r ^ tq the fluctuations are 
relatively weak, ya.Tu{\ip^{r)\) <^ {\i^^{^)\)u- In the opposite limit, r ^ tq, the fluctuations 
are expected to have essentially the GUE statistics, similarly to what we have found for the 
quasi-lD sample geometry, see eqs.( p.72D -( p.74|) . 



V. ASYMPTOTIC BEHAVIOR OF DISTRIBUTIONS OF VARIOUS 
QUANTITIES, AND THEIR INTERRELATIONSHIP. 

In the preceding sections, we have considered the spatial structure of ALS responsible 
for the asymptotic behavior of the distribution of eigenfunction intensity, V{u). One can 
consider, however, other quantities, large values of which indicate in some sense a stronger 
localization of an eigenfunction, so that the asymptotic behavior of corresponding distribu- 
tion function is also related to a kind of ALS. An illustrative example of such a quantity is 
the inverse partcipation ratio (IPR) I2 = ! d'^r\%l}'^{r)\. Let us consider the case of an infi- 
nite quasi-one- dimensional sample, for which the distribution function V{l2) was calculated 



exactly in Ref. [^,|T^]. It has the following asymptotic behavior: 

P(/2)~exp|-ye5/2| , (5.1) 

where we omitted preexponential factors. Let us suppose that the asymptotic form of V{l2) is 
determined by the same anomalously localized sates, which control the asymptotic behavior 
of V{y\). The intensity of such an ALS was found in Sec.|T| to be 

\^\r)\ = |^'(r)|g^oothl'^'(^)l ' « ^ (5-2) 
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where 

l^'(0lsmooth = ^(7^^ ' ief = 4^^ (5-3) 

and |$^(r)| exhibits the GUE-hke fluctuations ( p.74|) . The distribution V{u) behaves asymp- 
totically as V{u) ~ exp{— 4-\/uS^}; here the factor exp{— 2a/uS^} is the GUE probability of 
the quasi-jump in the vicinity of r = with 1^*^(0)1 = 2^JuS^. The remaining exp{— 
factor is therefore the weight of the envelope configuration IV^^I'^) Igniooth' '^^^ corresponding 



IPR is equal to I2 = {l/?>)\Ju/ S^. Thus, assuming that this ALS determines the asymptotics 
of the IPR distribution V{l2), we would get 

P(/2) ~exp {-6^5/2}. (5.4) 



This result has the same exponential form, as the correct asymptotics (|5.1|) , but with a larger 
numerical coefficient in the exponent. The explanation for this discrepancy is the following: 
an ALS, which is optimal for maximizing the local amplitude u = \ip'^{0)\, does not optimize 
the IPR. A detailed study, which will be published elsewhere |^^, shows that the ALS, 
which determine the asymptotics of the IPR distribution, have the following spatial shape: 

l^^(-)lsmooth = ^^Jf^^ ' = ^ ' (5-5) 



that is different from eq.( ^.3| ). Thus, the spatial form of an ALS depends on the specific 
physical quantity (local amplitude, IPR, . . . ), for which it represents an optimal fluctuation. 
In the remaining part of this section, we consider from this point of view the distributions 
of LDOS and of the relaxation times in open metallic samples, which have been studied via 



the supersymmetry approach in Refs. |22,23| and Ref. |24], respectively. 



A. Quasi-ID geometry. 

1. Distribution of relaxation times. 



The long-time dispersion of the average conductance has the LN form |^| 

G'(t)~exp|-f7ln2^-^| ; t»A-i, (5.6) 
where A = I/uqLS is the mean level spacing and g = 2tcvqDS/ L is the dimensionless 



conductance. We represent eq. (|5.6| ) as a superposition of the simple relaxation processes 



with mesoscopically distributed relaxation times [0]: 

Git) ~ j dt^e~"'^V{t^) (5.7) 
The distribution function V{t^) then behaves as follows: 

Vit^) ~ exp{-(7 \Ti\g/\t^)} ■ t't'^J^ (5-8) 
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This can be easily checked by substituting eg. {f).8\) into eq. (|5.7| ) and calculating the integral 
via the sationary point method; the stationary point equation being 



2gt^\n{gAt^) = t (5.9) 

Note that = gA is the inverse characteristic time of diffusion through the sample 
(Thouless energy), or in other words, the typical width of a level of an open system. The for- 
mula ( |5.8| ) concerns therefore the states with anomalously small widths in the energy space. 
The corresponding saddle-point configuration is found from the requirment of providing a 
minimum to the action 

s = I [rvstr(VQ)^ ^^_WD f \\v{yef , (5.10) 



A 

with the additional restriction 

t 



f d'^r [cosh e -1) = , (5.11) 

and the boundary condition on the boundary with leads ^'lig^^jg = 0. For a quasi-lD sample 
of the length L, the solution can be approximated as p2 



where 60 satisfies the equation 



= -tAOo , (5.13) 

TT 



SO that 

9o ^ In (^^tAln(iA)j ~ In (^^gAt^\n\gAt^j^ . (5.14) 

Relying on our previous experience, we believe that the smoothed intensity of the corre- 
sponding state is |'?/'^('^)lsixiooth stated in Ref. ||2^. Normalizing it by the 
condition J \ iP'^{t) \ = 1, we get 

I^^WIsmooth = ||e--oMA (5.15) 
Thus, the ALS, which gives a minimum to the level width t^^, has an exponential shape 



2. Distribution of local density of states. 

Now, let us consider the distribution V{p) of LDOS. Typically, in an open metallic sample 
the LDOS p{E,r) is given by a superposition of ~ l/toA = g adjacent levels, since their 
widths are of order of l/tz). However, we can expect that for p much greater than its average 
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value z/q, the asymptotic form of V{p) is determined by a probability to have a single narrow 
resonance, which gives this value of LDOS p{E,r). The most favorable situation happens 
when the resonance is located around the point r in the real space and around the energy 
E in energy space. The LDOS provided by such a resonance is expected to be: 

PALS = \^\r)\^ , (5.16) 

TT 

where t^^ is the resonance width. Thus, the optimal fluctuation shoul provide now a maxi- 
mum to the product of the local amplitude u = |^/'^(t*)| and the inverse level width t^. Since 
the distribution V{t^), eq. (|578|) , decays much slower than V{u), eq. (|3.36|) , one should expect 
the asymptotic behavior of V{p) to be mainly determined by P(t</,). Indeed, it was found in 
Ref. that V{p) has a LN form, similar to that of V{t^): 

Vip) ~ exp (^i- + -^^ ln2(p/z/o)} , (5.17) 

where, as before, L+ and L_ are the distances from the observation point r = to the 
sample edges. The corresponding saddle-point configuration reads: 

" 1 (pM)'-!'-!/^- , r<0 ^^-^^^ 

If we put the observation point in the middle of the sample, L_|_ = L_ = L/2, the configu- 
ration ( |5.18| ) acquires the same form as the optimal configuration ( p.l2| ) for the relaxation 
time t^. The corresponding values of and p are related as follows: 

-gAt^\n\gAt^) = p/uo , (5.19) 

TT 



or, expressing through p, 



(5.20) 



AgAuoln'ip/uo) 



Now, we calculate the value of the local amplitude |V'^(0)| for an ALS corresponding to the 
cofiguration (^.18|) . First, its smoothed intensity is given by 



^■■^('•)lsmooth-A^-'^""-!^(^)" (5.21) 

Second, the quasi-jump induced by the GUE-type fluctuations gives an additional factor, 
which can be found in the same way as prescribed by eq.( p.30 ): 



A{p) = -P^ \nV{p) = e} Hp/u,) (5.22) 

Combining eqs.( p.2(]| ), (|5.21|) , and ( ^.22|) , we can compute the LDOS (|5.16|) determined by 
this resonance state: 
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PALsiE^O) = |V^^(0)|3^,,th ■ Mp) ■ ^ = ■ 2gHp/.o) ■ ^^^^^^ = P (5-23) 

We have explicitly checked therefore that the LDOS p is indeed determined by a single 
ALS, smoothed intensity of which is given by eq. (|5.21|) . There are three sources of the 
enhancement of LDOS: i) amplitude of the smooth envelope of the wave function, ii) the 
short-scale GUE "bump", and iii) the inverse resonance width. They are represented by 
the three factors in eg. ( p.23| ), respectively. Note that the calculated pals reproduces the 
value of p with an amazing accuracy (including logarithmic factors and even the numerical 
coefficient). 

3. Distribution of global density of states. 

Finally, we discuss the contribution of ALS to the asymptotical behavior of the distribu- 
tion function 'P(i^) of the global density of states (DOS), 

u{E) = ^{Y.S{E - E^)) = i| d'rp{E,r) (5.24) 
A resonance state with an energy E and width t~^^ gives a following contribution to ^{E): 

^ALs{E) = -% = -UAuo (5.25) 

Thus, if we assume that the asymptotic behavior of 'P(i^) is determined by isolated (in energy 
space) anomalously localized states, it will have the form: 

V{u) ^V{t^ = ^) ~ exp [-g\n\gv/vo)} . (5.26) 

We will see below that an analogous procedure in 2D leads to a result for V{v) which is in 
full agreement with the renormalization group calculation of Altshuler, Kravtsov and Lerner 



B. 2D geometry. 

1. Distribution of relaxation times. 



Distribution of relaxation times has the form 23 



with g = 2ttvqD being the dimensionless conductance of a 2D square. An ALS corresponding 
to the first regime, t£, <^t^ <^tD{L/l)^, has the following spatial structure: 



26 



IQnDt^ [(r/L)2 + Ly{16Dt^)]^ 



(5.28) 



so that it has an effective locahzation length ^e/ ~ L(t£)/t^Y^'^, with the intensity decreasing 
as l/r"^ outside the region of the extent C,ef- As to the ultra-long-time region, tn^L/l)"^ , 
the saddle-point solution reads: 



od{r) 



{r/L) 



7t-2 



where 



l^^ = Itl ; It 



2 ' 



<r<L 



(5.29) 



ln(t0/tg) 
ln(i?/0 



(5.30) 



Now, It plays a role of an effective localization length, and the intensity shows the following 
behavior: 



mr)\ 



r 



-7t-2 



7W2 I ,(*) 



<r<L. 



(5.31) 



2. Distribution of local density of states. 
Distribution of LDOS, V{p), has the following asymptotics 



V{p) ~ exp 



tt^z/qD In p ' 



The corresponding saddle-point solution reads 



uo \ r 



where li^^ = jpl, and 



7p 



ln(p/z/o) 



Normalizing it, we get the following ALS intensity at r > li''^: 



I ^^('^) I smooth 



(2 - 7p) 



(p)~ 



l_p_ 
V uo \ r 



7p 



7p<2 
7p>2 



(5.32) 



(5.33) 



(5.34) 



The value of the GUE-type quasi-jump is 
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A{p) = -p^\nV{p) = e^pn^iy,V (5.35) 

To estimate the escape time from this resonance state, we note that its power-law decay, 
[■^^(r)! oc {r /li^^)~'^p , is similar to that of the ALS optimizing the relaxation time, |'?/'^(r)| cx 
[r . This allows us to identify 7p = 7t + 2, so that 

mtD ~ <; ^ (z) ' ^ ^ (5.36) 

1 , 7p < 2 , 



up to logarithmic prefactors depending on 7p and 7^. Substituting eqs.( ^.34D , ( ^.35|) , and 



p.3(j| ) in eq.( |5.16| ), we get Pals ~ P- Thus, we have checked that the ALS determined by the 



saddle-point solution (|5.33|) indeed provides the value of LDOS which is equal (within the 
accuracy of our consideration) to p. This confirms an assumption that an anomalously large 
value of p is typically governed by a single ALS with the spatial structure described by the 
corresponding saddle point configuration. Like in the quasi-lD case, the enhancement of p 
is determined by the product of three factors: |V'^(^)lsmooth' ^(p)' t^p/tD, represented 
by eqs. (|5.34|) , (|5.35| ), and (|5.36|) , respectively. 



3. Distribution of global density of states. 

Now, we consider the contribution of ALS to the asymptotics of the distribution of global 
DOS, which can be estimated according to eq. (|5.26|) as follows: 



^M~^('* = ^)~],,J ^g lnVMAr) ! (5.37) 

I 2 ln(L//) J 



I/O g \ l 



The far LN asymptotic tail in eq. ( |5.37| ) is in full agreement with the RG calculation by 



Altshuler, Kravtsov, and Lerner 0. We find also an intermediate power-law behavior, 
which could not be obtained from the study of cumulants in Ref. [0. We note, however, 
that this power-law form is fully consistent with the change of the behavior of cumulants 
at n ~ TT^f discovered in 0. 
Finally, taking into account the close similarity between the cumulants of DOS and of 
conductance 0, it is natural to suppose that the distribution function of conductance in a 
2D metallic system has essentially the same behavior ( |5.37|) , with an intermediate power- 



law regime. This hypothesis, which would be fully consistent with a power-law asymptotic 
behavior of conductance distribution function on the mobility edge suggested by Shapiro 



34 , needs however further verification. 



C. 3D geometry. 



As has been already mentioned in Sec.|rV|, in three dimensions the ALS are in fact not 
localized. Their intensity just shows a relatively "narrow" bump on top of the usual average 
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value IV'^I = Comparing the results of Sec.|^ and of Refs. [^,|^, we find that the 

spatial shape of these "bumps" is the same for the states, which are optimal for all the 
distributions V{u), V{p), and G{t). Namely, it has the form 



l^^(^)lsmooth ~ ^^^P ' i = u,p,t, (5.38) 

where 

uV 
t 

and Ci are numerical constants of order of unity. The corresponding distributions have one 
and the same, log-cube-exponential, form: 

V{i) ~ exp{-const(A;i7/)^ In^ Zi] (5.39) 

Note that formation of a large value of LDOS p{E,r) cannot be explained in 3D as a 
contribution of a single ALS. Indeed, in a metallic sample LDOS is typically provided by a 
number of levels of order of g ^ kplL. In order that a single level might give such (or even 
larger) a value of LDOS, it should has a local amplitude (or, alternatively, an inverse width 
in energy space) enhanced by a factor of g. However, in contrast to the quasi-lD and 2D 
situations, g does not enter the asymptotics ( |5.39|) , which do not depend on the system size 
L. Therefore, a high value of LDOS is in 3D typically due to contribution of a large number 
(oc L/l) of adjacent levels. 



Zu — 
Zt = 



VI. ADDITIONAL COMMENTS. 

A. States localized near the boundary. 

We assumed throughout the paper that the center of an ALS is located far enough from 
the sample edge. For a quasi-lD sample, this means that C,ef Lj^, L_. In the 2D case this 
implies that the distance from the observation point to the boundary is of the same order 
of magnitude in all directions, so that ln(L//) is defined without ambiguity. Here, we will 
consider briefly the role of ALS situated close to the boundary, when these conditions are 
violated. 

We start from the quasi-lD geometry. Let us calculate the distribution function V{u) in 
a point located very close to one of the sample edges. Formally, this means that L_ <C ^e/- 
Then the function W^^\uS^, r_) in eq.( |3.4| ) can be approximated by unity, and we get 



V{u) = le^'S'/'L-'/\--'/^exp I -2^+^ ( 1 - + ...\\ (6.1) 
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We see therefore that close to the boundary the distribution V{n) has the asymptotic 
decay V{u) ~ exp{—2^/uS^}, which is slower than in the bulk of the sample, V{u) ~ 
exp{—A^/uS^} . This means that if we consider the distribution V{u) averaged over the po- 
sition of the observation point, its asymptotic tail will be always dominated by contribution 
of the points located close to the boundary, V{u) ~ exp{— 2a/uS^}. This could be already 
anticipated from eg. ( ^.461 ), where the factor exp |^ (^-^ + 77-)} strongly increase with ap- 
proaching one of the sample edges. The same tendency, but in a weaker form, is observed in 
eqs.( p.lOD , ( |2.11| ), ( p.l2| ). Calculating the average intensity {\4''^{r)\)u of the corresponding 



ALS, we find that at r > / eqs.( |3.5(]|) , ( p.54| ) retain their validity, with an additional overall 
factor of 2. At small r, eq.( p.44| ) is slightly modified: 



\ 1/2 



(6.2) 



In 2D, we can consider a sample of the semicircular shape, with the observation point 
located in the center of the diameter serving as a boundary. The saddle-point solution then 
has exactly the same form (f4.2|) , and the ALS intensity is still given by eq. ( |4.5|) , with an 



additional factor 2. The asymptotic form of the distribution function V{u) gets an extra 
factor 1/2 in the exponent: 

We expect this result to be applicable to any 2D sample of a characteristic size L, with a 
smooth boundary and the observation point taken in the vicinity of the boundary. 

We see therefore, that, very generally, the probability of formation of an ALS with the 
center in a given point is strongly enhanced (via an extra factor 1/2 in the exponent), if this 
point lies close to the sample edge. This leads to the additional factor 1/2 in the exponent 
in the asymptotical form of the distributions V{u) and V{p) near the boundary. 



B. Orthogonal symmetry class. 

All the considerations in this paper can be straightforwardly generalized to the systems 
with unbroken time reversal invariance (orthogonal symmetry class). The main results are 
as follows: 

i) all formulas for the average spatial density {\'ip'^{r)\) in the metallic samples retain their 
validity. In particular, in the quasi-lD case, eqs.( p.43D , ( p.44| ), (|3.5CI| ), (|3.54|) hold with 
the same definition of ^ = 211^^30. In the far localized tail, eq.( |3.42| ), ^ is replaced 
by (^/2, which is just the conventional dependence of the localization length on the 
symmetry of ensemble; 

ii) in the expressions for the asymptotics of all distribution functions, an extra factor 1/2 
appears in the exponent; 

iii) GUE-type fluctuations are replaced by the GOE-type ones, where appropriate. 
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VII. SUMMARY. 



In this paper, we have studied the spatial structure of the anomalously localized sates in 
weakly disordered samples. Such states appear to govern the asymptotical behavior of the 
distribution function V{u) of local amplitudes of eigenfunctions. In the quasi- ID geometry, 
an ALS has an effective localization length C,ef much shorter than the conventional one ^. 
We were able to calculate exactly the average intensity of such a state. The found spatial 
distribution of the ALS intensity, {\4''^{r)\) oc l/(r + ^e/Y, turned out to be in agreement 



with the form of the solution of the saddle-point equation of Ref . |21| . Thus, the saddle-point 



configuration indeed describes the average intensity of an ALS, as was conjectured in Refs. 



2^ , |2T|] . This allowed us to describe the spatial structure of ALS in 2D and 3D, as well. We 



have also studied the fluctuations of the ALS intensity and found them to be essentially of 
the same GUE type, as for a typical delocalized eigenstate. 

Not only the asymptotic behavior of V{u), but also that of distributions of other quan- 
tities can be governed by a kind of ALS. Here are several important examples of such 
quantitites: inverse participation ratio I2, relaxation time t^j,, local density of states p{E,r), 
global density of states i^{E). We have found that the spatial structure of ALS relevant 
to the asymptotic behavior of different distributions may be different. This is because an 
ALS constitutes an optimal fluctuation for one of the above quantities, and the form of 
this fluctuation depends on the specific characteristic, which is to be optimized. Finally, we 
have discussed interrelations between asymptotics of various distributions mentioned above. 
In the quasi-lD and 2D cases, this allowed us to present a comprehensive picture, which 
explains all the asymptotics as governed by exponentially rare events of formation of ALS. 
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APPENDIX. JOINT DISTRIBUTION FUNCTION OF EIGENFUNCTION 
INTENSITIES IN TWO DIFFERENT SPATIAL POINTS. 

In this Appendix, we write down the results for the joint distribution function of inten- 
sities of an eigenfunction in two spatial points, 

p(n, V) = {6mr\ - n)smvf\ - p) , (a.i) 

for a quasi-lD system. For / < r <C ^, the function V{u, v) can be restored from its moments, 
eq. (|3.61 ). The result is obtained in the following form: 



V{u,v-r) = J--^^ rdi/z/sinh(7rz/)e--i(i+'^')/^ f' d^^^W^'^ (-^,^2 



31 



In the opposite case, r < I, the expression ( p.75| ) for the moments has to be used, yielding 
z{u, V, r, (f)) = ^S \ + 



1 + JkJr)e^^ 1 + JkJr)e-^'t> 



Unfortunately, even in the simplest case of an infinitely long sample, when the function W^^"* 
is given by eq. (|3.32|) , formulas ( [A .21) and ( |A.3D are too involved. This is why we have chosen 
in Sec.|TT| to analyze the expressions for the moments, rather than the distribution function 
itself. 
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